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Abstract 

This article gives a pedagogic derivation of the Bethe Ansatz solution for ID interacting anyons. 
This includes a demonstration of the subtle role of the anyonic phases in the Bethe Ansatz arising 
from the anyonic commutation relations. The thermodynamic Bethe Ansatz equations defining 
the temperature dependent properties of the model are also derived, from which some groundstate 
properties are obtained. 



PACS numbers: 02.30.Ik, 05.30.Pr 



I. INTRODUCTION 



A number of ID models in quantum many-body physics have been solved by the Bethe 
Ansatz following Bethe's pioneering work on the exact solution of the ID Heisenberg mag- 
netic spin chain in 1931 . During the past 75 years the Bethe Ansatz has been developed 
and applied to physical problems such as the ID (5-function interacting Bose j^, 4] and Fermi 
jsl gases, the ID Hubbard model and 2D vertex models 7, 8| (many of these key papers 
are collected in Ref. ^]). The solution of such models contributed to the development of the 
celebrated Yang-Baxter equation, which gives the consistency conditions for many-body scat- 



tering problems and plays a crucial role in quan tum integrable systems 



and in modern mathematical 
inverse scattering method 11 
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pnysics Il7l |. The related establishment of the quantum 



llSl . Il9l | provided widespread applications of t 



le Yang-Baxter 



equation in low-dimensional quantum systems, such as the Kondo problem 



son model 



20|, the Ander- 



2l| and long range interaction systems 22l. |23|. 



A further pioneering application of the Bethe Ansatz was to the ID (5-function interacting 
Bose gas [sl which is related to the quantum nonlinear Schrodinger equation. This model 
provides an important realistic physical description of an interacting ID Bose gas. It is 
also arguably one of the most simplest and pedagogic models solved in terms of the Bethe 
Ansatz. In general the applicability of the Bethe Ansatz depends on the reducibility of the 
multi-particle scattering matrix to the product of many two-particle scattering matrices. The 
starting point for the Bethe Ansatz approach in quantum many-body physics is to reduce the 
eigenvalue problem of the field theoretic hamiltonian into a quantum mechanical many-body 
problem. The wavefunction of the many-body hamiltonian inherits the statistical signature 
of the interacting particles, which leads to striking and subtle quantum many-body effects. 
For example, significantly different quantum effects between the ID 5-function interacting 
Bose and Fermi gases are seen clearly from the Bethe Ansatz solutions 0, Isj. 

On the other hand, anyons 2J, |25|] may also exist in both two and one dimension, obeying 
fractional statistics. For a 2D electron gas in the fractional quantum Hall (FQH) regime, the 
quasi-particles are charged anyons A more general description of quantum statistics is 
provided by Haldane exclusion statistics {2^], which is a formulation of fractional statistics 
3ased on a generalized Pauli exclusion principle, now called generalized exclusion statistics 
28l . |29| . In ID, a wavefunction with anyonic symmetry may also occur [30l]. Anyons in ID 
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acquire a multi-step function-like phase when two identical particles exchange their positions 



in the scattering process. This topologica 



interacting model of anyons 
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34 



3hase results in rich quantum effects in the ID 



35|. 



In this paper, celebrating the 75th anniversary of the Bethe Ansatz, we give a detailed 
and pedagogic derivation of the Bethe Ansatz solution for ID (5-function interacting anyons, 
discussed for the first time by Kundu [scl- The ground state properties derived from the 
thermodynamic Bethe Ansatz are also presented. 

This paper is set out as follows. In section [Tll we introduce the ID interacting anyon 
model and show the reducibility of the field-theoretic Hamiltonian to a quantum mechanics 
many-body problem. The Bethe Ansatz solution is derived in section IIIII In section \1V\ we 
study the thermodynamic Bethe Ansatz equations (TBA). Section|V]is devoted to concluding 
remarks. 



II. THE MODEL 

We consider N anyons with 5-function interaction in ID with hamiltonian [3^ 

fi^ 1 

H = / dxd'^Ux)d'^(x) + - giB / dx "^UxW (xW (xW (x) (1) 

2m Jq 2 Jo 

under periodic boundary conditions, with x a coordinate in length L. Here m denotes the 
atomic mass. Hereafter we set h = 2m = 1 for convenience. The coupling constant qib 
is determined by giD = ti^c/m where, at least for the Bose gas, the coupling strength c is 
tuned through an effective ID scattering length air, via confinement in experiments. We 
also use a dimensionless coupling constant 7 = c/n to characterize different physical regimes 
of the anyon gas, where n = N/L is the linear density. The operators and \l/(a:) 

are the creation and annihilation operators (or operator valued particle density) at point x 
satisfying the anyonic (equal-time) commutation relations 

^{xi)^^{x2) = e-^'^^(^i'^2)$t(a;2)^(xi) + 5(xi -X2) 
^(xi)^(x2) = e'''"'("i'"2)^(x2)^(xi), 
^^ixi)¥{x2) = e''^"'("i'"2)^^(x2)^^(xi). (2) 

Here the multi-step function w(xi, X2) = —w{x2, xi) = 1 for order xi > X2, with w{x, x) = 0. 
The anyonic phase w{x, x) = for two colliding particles has bosonic signature at the point 
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Xi = X2- We note that in imposing the restriction e^'^'^^^'^) = e''^ on the anyonic phase the 
above commutation relations give hard-core relations which have been constructed based on 
lattice integrable models 36| and also suggested in the anyon- Fermi mapping for the anyonic 



Tonks-Girardeau gas 3J]. We shall see that this phase is not suitable for the continuous 
model of anyons. In contrast to the ID interacting Bose gas js], hamiltonian ([1]) exhibits 
both anyonic statistical and dynamical interactions, which can map to a ID interacting Bose 
gas with multi-(5-function and momentum- dependent interactions 3^. The more general 
anyonic statistical and dynamical interactions result in a richer range of quantum effects 
than those of the ID interacting Bose gas. 

We first consider the corresponding equation of motion —idt'^{x^t) = [H, \E'(x, t)] via the 
time-dependent quantum fields, i.e., the nonlinear Schrodinger equation 



(3) 



In order to appreciate the subtlety of the anyonic phase, the nonlinear Schrodinger equation 
can be calculated as 

£dx^ (^^,,^\x^)^,,^{x^) + ^g^J,{^^ {xi))\^{xi)f^ ,^(x2) 
= £ dx, (^-^\xi)di^^{xi) + ^(?iD^^'(xi)^'(a;i)^ , ^(0:2) 

= - [ dxi (e''^"'(^i'^2)^t(xi)^(x2) - e-'^'"^^^^^'^^^{xi)'^{x2) - S{x2 - Xi)) dl^^{xi) 
Jo 

+ IgiD j dxi [^^\x^)^\xi)^{x2) - {e-''^^^''''''^'^^^\xi)^{x2) 

+ e-'"'"("^'"^)5(x2 - xi)¥{xi) + e-=''""'("^'"i)5(x2 - xi)^\xi)) ^^{xi)] 



dl^{x2)-giD'^\x2)^\x2). 



(4) 



In the above equation the time t has been omitted for simplicity. The properties w{x, x) = 
and w{xi,X2) = —w{x2, Xi) have been used. We stress that the phases induced from exchang- 
ing operator positions in the terms '^{xi)'^{x2) and '^{x2)'^^{xi) are some fixed constants, 
depending on the relation between xi and X2. The derivative operator d^^ then only acts 
on the field operators. Furthermore, if we choose e^'^'"^^'^) = e"^, the last term in the above 
equation does not exhist due to the cancellation of the terms e~^'^'^^^^'^'^H{x2 — a;i)\I'^(xi) 
and e-2i«"'(^2,:!:i)^(a,2 - xi)'^^{xi). This means that there is no 5-function interaction in the 
quantum mechanical many-body hamiltonian without internal spin degrees of freedom due 



to the hard-core behaviour. For this reason we take the choice of phase w{x, x) = in order 
to ensure the commutation relations ([2]). Under this requirement we can obtain the form of 
the nonhnear Schrodinger equation ([3]). 



Following the way suggested for the ID Bose gas 



we define a Fock vacuum state 



\E'(x) |0) = 0. Thus the number operator N and momentum operator P are 



L 



N = / «^(a;)^(a;), (5) 







P = i /" dx [d^'^^ix)] ^(x). (6) 
Jo 

To properly denote the anyonic phases in the eigenfunctions | $) of hamiltonian ([1]), we 
assign particle coordinates Xi in the order xi < X2 < ■ ■ ■ < xn. Based on this assigned order 
the multi-step function w{xi, xj) can be used to count phase alternation when two particles 
are interchanged, i.e., w{xi,Xj) = —w{xj,Xi). The iV-particle eigenstate is written as 

|$)= / dx''e-''^x{xi...XN)'^\xi)...^\xN)\0) (7) 
Jo 

where the Bethe Ansatz wavefunction is of the form 

Xixi ...Xn) = e"^ ^^»<-^ '"^'""'''^ ^i^P^ ■ ■ ■ A;p7v)e^^*^"^i+-+'=^^^^). (8) 

p 

Here the sum extends over all A^! permutations P. The order in which the particles are 
created in the A^-particle eigenstate incurs the phase factor in the wave function (IHl). 
Integration involves changes in the order of creating particles due to the permutation of 
coordinates. We can denote the permutation P as 



P 



1 2 ... A^ 
PI P2 ... PN ^ 

which transforms the m-th object into the Pm-th position for 1 < m < A^ In general, 
for two different permutations we have 








QP 



-1 



PI P2 ... PN 

Qi Q2 ■ ■ ■ Qn 



In order to understand the above expressions, we give the exphcit form of the eigenstate 
([7]) for = 3 in Appendix [Aj The coefficients A{kpi . . . kp^) are obtained exphcitly via 
the Bethe Ansatz approach in the next section. By comparing eqns (lA2p and (lA4p we see 
clearly that the phase factors in the multi-valued wavefunction ([HD are diminished by those 
from permutations of the particles in the eigenstate 1$) such that the integrand in ([7]) is 
single valued and fully symmetric. We extracted a global phase factor q-^'^^/'^ in order to 
symmetrize the anyonic phase factor in the wave function ([8]) so that it has k ^ k + 47r 
symmetry. The eigenstate still has k ^ k + 27r symmetry. Moreover, it is easily seen that 
the wave function satisfies the anyonic symmetry jsO] 

x{...Xi...Xj...)= e-'Kn=.+i^(^-^'=)-n=-+i"'(^.-^'=))^(. ..Xj...x,...). (9) 

Acting on the eigenstate ([7]) with the operator P defined in (E]) gives 

P I $) = i / dy r dx^e-'^^xixi ...xn) 
Jo Jo 

N 

X J2 e-«Sti^ -(J^-->)vI/t(xi) . . . vl>t(x,_i)5(z/ - Xk)^^{xk+i) . . . vl>t(x^) | 0) 

k=l 
j-L 

= i / dx^e~^~x{xi ■ ■ ■ ^n) 
Jo 

N 

X ^^(a^i) • • • '^\xk-i)d,^^\xk)^\xk+i) . . . ^H^n) I 0) 

k=l 

= £ dx^'e-'^ I -i ^ d^.xixi ...xn)^ ^\xi) . . . ^\xn) 1 0) 



(10) 



It therefore follows that we can obtain the quantum mechanical momentum operator as 

N 



iE^^- (11) 



-1 

k=l 



Similarly, by acting on the eigenstate ([7]) with the hamiltonian ([T]) the eigenvalue problem 
for hamiltonian ([1]), namely 

H\^)= dx^e-'^HNx{xi...XN)^\xi)...^\xN)\Q), (12) 
Jo 

can be reduced to solving the quantum mechanical problem 

Hnx{xi- ■ -Xn) = Ex{xi. . .xn), (13) 
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where 



^ Q2 

^^ = -^J2d^+3iD Yl ^i^i-^j)- (14) 

i=l * i<«<i<JV 

The quantum mechanical hamiltonian (fT^ describes the ID (5-function interacting quantum 
gas of anyons confined in a periodic length L. The details of the above calculation are 
lengthy but straightforward. Without loss of generality, we give the details of the proof for 
the case = 3 in Appendix [Bl Again, the subtlety of the anyonic phases is such that they 
just match the applicability of the Bethe Ansatz - from the eigenvalue problem of the field 
theoretic hamiltonian ([1]) to the quantum many-body problem (fT^ . The next task is to 
determine the wavefunction by means of the Bethe Ansatz. 

III. BETHE ANSATZ SOLUTION 

In the preceding section, we showed the equivalence between the quantum field theoretic 
and quantum many-body problems. The Bethe Ansatz wavefunction ([8]) was written for 



32l |. The wavefunction in 



the fundamental domain < Xi < X2 < . . . < xjy < L 
other domains is an extension of the wavefunction ([H]) via the anyonic symmetry (Q. The 
(5-function potential causes a jump in the derivative of the wavefunction in the eigenvalue 
equation f|T3|) . Changing the coordinates to centre of mass coordinates X = {xj + Xk)/^ and 
Y = Xj — Xk leads to the eigenvalue equation (fT3|) in the form 



Integrating both sides of this equation with respect to Y from — e to +e and taking the limit 
e — gives the discontinuity condition 



(f^Xj ^Xi) \{X\) • • • ; Xi^ Xj, . . . , Xj\f^ \xj=Xi+e ij^Xj ^Xi) X.{X\, . . . , Xj, Xi, ■ ■ ■ , Xj\f^ |a 
= 2cx{Xi, . . .,Xi,Xj, . . .,XN)\xi=xy (16) 

on the derivative of the wavefunction. The two particle scattering relation can then be 
worked out directly from the wavefunction xi^i^ . . . ,xn)- This model does not have internal 
degrees of freedom. Therefore, it is not necessary to introduce extra coordinate indices to 
distinguish the coefficients in the wave function for different domains. We may write the 
wavefunction as 



X{xi ...xn) = e-'^("("-"^-)+^^<-'" [. . . A(. . . ki, kj . . .)e'(- 



■ Xj -|-. . .) 



+A{. . . kj, ki . . ^ _ _ _j 



(17) 



for the domain < Xi < X2 < ■ ■ ■ Xi < Xj . . . < < L and 



X{xi ...xn) 



N 



^("''"'")) [...A{...ki,kj...)e"^- 



■ ~\~k'iXj-\-hj -|-. . .) 



+...] 



(18) 



for the domain < Xi < X2 < ■ ■ ■ Xj < Xi . . . < x^ < L. For Xj = Xj, the wavefunction is 
given by 



-|-y4.(. . . kj, ki . . .)e 



i{... + {ki+kj)xi+...) 



+ ...]. 



(19) 



Apphcation of the derivative discontinuity condition (|T6l) on the wavefunctions f|T71) . f|T8|) 
and (fT9|) gives 

■.2c[A{...ki,k,...) + A{...kj,ki...)]. (20) 



This equation gives a relation between the coefficients A{kpi . . . kpjq) of the form [30|, l32| 

~\~ ic 



kj kj ic 



• • ki^ kj . . .) 



(21) 



which is the two-body scattering relation. Here the anyonic parameter k and the dynamical 
interaction c are inextricably related via the effective coupling constant 30|, l32| 



c' = c/ cos(k/2). 



(22) 



This results in an interesting resonance-like effect in the effective coupling constant c' with 
respect to the statistical interaction around k = ir 3^. The single- valued eigenstate ([7]) 
makes the continuity of the wavefunction ([8]) well defined because the phase factors in the 
multi-valued wavefunction ([8]) can be cancelled by those from permutations of the particles 
in the eigenstate | $) for different domains. The probability density can be written as 
|x(a;i, . . . , xat)!^ which is subject to the scattering problem. It suffices to request a continuity 
relation of the wavefunction of the form 



lim xi^i, ■ ■ ■,Xi,Xj, ...,xn) 



\im x{xi, ...,Xj,Xi,..., xn) = x(xi, . . . ,Xi,Xi, . . . ,xn), 



(23) 



or equivalently 

lim [A{. . . ki, kj . . _y(-+k.x,+kjx,+...) + _ _ ^.^ ^. _ _ _)gi(...+fc,x,+fc,a;,+...) ^ _ _ _j 

Xi — Xj 

= lim [A{. . . ki,kj . . .)e'^---+''^''^+''^'''+---^ + A{. . . kj,ki . . .)e'^^^^^^^ (24) 

which leads to an identity equation for the coefficients in the above two domains. 

The periodic boundary conditions xi^i = 0,X2 ■ ■ ■ xn) = x{^2 ■ ■ ■ xn,xi = L) for the 
system are equivalent to x{^i^^2 ■ ■ ■ xn) = e~^'^^^~^\{xi + L, X2 ■ ■ ■ x^)- The anyonic phases 
for Xi at the left end and at the right end of the domain are different due to the changes in 
the order of particles. Explicitly, 

Q-f{Y.Z2^i^i'''i)+T,i<i<jw{xi,xj)) ^ A{kpi ■ ■ ■ kpN)e^^^^^'^^""^'"'^'''^^ 

p 

Q 

Notice that the second sum in the phase factors in both sides will cancel in deriving the 
Bethe ansatz equations below. In order to compare the coefficients on both sides of the 
above equation, we permute the order Q back to the order P. There are A^! patterns related 
to all different ways of displaying the positions of fc's. For example, in the right hand side 
of equation fl25|) we can permute ki back from the right side of 

ki k2 ... ^Af-l kjy 

k2 ks . . . kN ki 

to the left side 

ki k2 ... kN-i kN 
ki k2 ... kN-i kN 
The Bethe Ansatz equations (BAE) 



i=i ■> ^ 

follow directly through the scattering process. We demonstrate the derivation of the BAE 
fl26|) in Appendix O for = 3. 



It should be noted that this model reduces to the interacting Bose gas 3" for k = 0. For 
K = 7T and Sir it reduces to the free Fermi gas. The anyonic gas lies in the range < /t < vr 
and 37r < K < An, where the effective interaction c' > 0. However, if the anyon parameter 
K is tuned smoothly from k, < tt to n > ir, i.e., vr < k < 3tt, the effective interaction is 
attractive. The ground state properties of this model have been studied in Ref. [3^ . 



IV. THERMODYNAMIC BETHE ANSATZ 



In general the exponential phase factor in the BAE (126!) . picked up from the statistical 
interaction during the scattering process, may shift the system into higher excited states. 
The total momentum is p = N{N — 1)k,/L + 2d7T/L, where d is an arbitrary integer. This 
can be easily seen from the logarithmic form 

^ k - k 

kjL = 27rlj + k{N -1)-2Y^ arctan( \ ' ) (27) 

1=1 ^ 

of the BAE, where Ij = j — {N+l)/2 with j = 1, . . . , A^. The energy is given hj E = Ylf=i ^'j- 
In minimizing the energy we consider k,{N — 1) = u (mod 2tt) in the phase factor with 
—TT < ly < 7T. Each quasimomentum kj shifts to kj + u/L in the ground state. In general the 
anyonic factor k, produces an effect like the self-sustained Aharonov-Bohm flux resulting in 
metastable states. In the thermodynamic limit, the BAE (|26|) can be written in the form of 
an integral equation 

Here Q is the cut-off momentum for the ground state. We have introduced the density of 
roots p{kj) = 1/ {L{kj^i — kj)) in the interval /Cj+i — kj. It follows that in the thermodynamic 
limit the lowest energy is given hj E = N{n'^e{'-f, k) + u'^/L'^) where 

e(7, k) = — J g{x)x^dx. (29) 

The last term in the energy E can be ignored compared to the kinetic energy. The root 
density g{x) := p{Qx) and the parameter \ = cjQ are determined by Lieb-Liniger type 
integral equations of the form 

1 Acos(k/2) g{.y)dy 

g[x) 



2tx vr 7_i A2 + cos2(|)(x - ?/)2 



A = 7 / g{x)dx. (30) 



The numerical solution of this equation has been given in Ref. [32 1. 

At zero temperature, the quantum numbers Ij form a uniform lattice from — (A^ — l)/2 
to (A^ — l)/2. The quasimomenta kj form a non-uniform distribution p{k) between k = —Q 
and k = Q, which is determined by the BAE (!28l) . The occupied lattice sites are called 
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quasiparticles. For the excited states the lattice sites are not all occupied. This means 
that some /c's move out of the Dirac sea leaving some lattice sites unoccupied, the so called 
holes. If one k leaves the Dirac sea, the remaining roots all move. This collective behaviour 
phenomenon is governed by the linear dispersion relation E — Eq ^ VcP as p — > 0, where the 
sound velocity is Vc ~ 2n7r(l — 47^-*^ cos(k/2)) [32i]. In this way, the density of quasiparticles 
and the density of the holes are determined by the equation 

1 1 f°° 2c'p(k')dk' 

The energy per particle for the state is given by 

rp poo 

— = n-^ p{k)k^dk (32) 
' J —00 

where the linear density n = J^^p{k)dk. 

To obtain the thermodynamics following the Yang- Yang approach 37(1, an equilibrium 

state associated with the entropy 

S 



n-' / [{p{k) + ph(A;)) Hp{k) + p^{k)) - p{k) In p{k) - p^{k) In p^{k)] dk (33) 
' J —00 

is a mixture of some eigenstates of the Hamiltonian. For the equilibrium state the density of 
quasiparticles is determined by minimizing the free energy F = E — TS — pN, i.e. 6F = 0, 
where p is the chemical potential. It follows that the TBA equation and the free energy are 
given by 

e{k) = e^{k)-p / d}^e\k-}^)\Yi{l + e-^—) (34) 

— = pn / cifcln(l + e"^) (35) 

L 271 

respectively. Here 

9'{x) = , 2ccos(./2) 

^ ^ C2 + COs2(k/2)x2 ^ ^ 

Here e{k) is the dressed energy e^^''^^'^ := p^{k)/ p{k) and the function e°(fc) = {k + u/LY. 



331]. For K = 0, the thermo 



38 



39] 



The thermodynamics of the system has been discussed in Ref. 
dynamics of the ID Bose gas has been studied extensively, see, e.g. Refs 

The TBA equation also provides ground state properties and a clear picture of the 
band filling. At zero temperature, the negative part e~ of the dressed energy makes a 
contribution to the free energy, i.e. for T — > 0, the TBA equation ( |34l) becomes 

11 



The pressure per unit length is given by Pq = I-q^ {k)dk. For the Tonks- Girardeau 
regime, i.e. for 7 ^ 1, we have 

'W^'^'-l^-^- (38) 

From this equation, with the help of the conditions e(Q) = and n = dPo/dfi, one can obtain 
physical quantities, such as the ground state energy Eq, chemical potential /i, pressure Pq 
and the cut-off momentum Q. The results are 

^ ,^ 1 3 2/ Acos{k/2)\ 2 2 A 16cos(k/2) 



3 \ 7 / \ 37 

202/ 6cos(/«/2)\ ^ / 2cos(/«/2)\ 
-o-g-Vj^l ^j, Q^nn[l l^j. (39) 



The macroscopic velocity is 



BP 

2-^ ^ 2mr{l - 47"^ cos(fi:/2)). (40) 



On the other hand, from the BAE (!28|) . we have 

where p is the total momentum. Then from the condition n = j3) p{k)dk, we have 



Q « ,„ ( 1 - !£2!i^/!)) , B„/i « i„v (1 - + ) (42) 



which coincide with the TBA results ( l39l) . 



V. CONCLUSION 

In this paper we have derived the Bethe Ansatz solution of the ID (5-function interacting 
model of anyons, which is a natural extension of the ID ^-function interacting Bose gas, 
recovered in the limit k = 0. The ground state properties fl39|) have been derived from the 
TBA equations fl34l) . The subtlety of the anyonic symmetry results in rich quantum effects 



for the degenerate anyon gas [32, l33| . Here we remark that a generalization of this model to 
the ID 5-function interacting model of anyons with internal spin degrees of freedom should 
be possible. We expect that the combination of anyonic statistical interaction, dynamical 
interaction and spin degrees will result in new quantum effects. 
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APPENDIX A: THREE PARTICLE EIGENSTATE 



Here we give the explicit eigenstate ([7]) for = 3. For the order of particle coordinates 
Xi < X2 < X3 the wavefunction is 



^-i'^{w{xi,X2)+w(xi,X3)+w{x2,X3)) 



J2 Mkpi, kp2, fcp3)e'=^i^^+'=^2^2+fcp3X3_ (Al) 



The sum runs over all 3! permutations P. According to the definition in the anyonic commu- 

. 3 

tation relations the above phase factor is e's'^ for = 3 in the assigned order xi < X2 < x^. 
In this domain the eigenstate for hamiltonian ([T]) is 



1$) 



^rf,^Q-^^Q-l§{w{xi,X2)+w{xi,X3)+w{x2,X3)) 

X^\XI)¥{X2)'^\X3)\0) 



+kp2X2+kp3X3 



dx 



J2 Mkpi, kp2, A;,3)e'="^i+'=-^^+'=^3^^ 



<i/\xi)¥{x2)'^\x3)\0). (A2) 



The phase factor is unity in the assigned order xi < X2 < x^. 

Now consider two-particle exchange, i.e. the particle order X2 < xi < X3. We can write 
the wavefunction as 



Xi^2 ^1 2^3) = Q-i'^(w(x2,xi)+w{x2,X3)+w{xi,X3)) 

Thus the eigenstate is given by 

X^^{X2)'^\XI)^^{X3)\0) 

^-i'^(w(xi,X2)+w{x2,X3)+w{xi,X3)) 



A{kpi, kp2, kp3)e 



kpiX2+kp2Xi+kp3X3 



(A3) 



J2 Mkpu kp2, fcp3)e''^^^^+^^^^^+''^^^^ 



J2 Mkpi, kp2, A;p3)e*^^i^2+'=^2^i+'^^^^^ 



dx^ 



J2 Mkpi, kp2, A;p3))e'=^i^2+'^^2Xi+fcp3X3 



^t(xi)^^(x2)^^(x3)|0). (A4) 



It is clear that the summation in the phase factor of the wavefunction x simply counts the 
changes of the order in the relevent domain. 
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APPENDIX B: FROM QUANTUM FIELD TO MANY-BODY QUANTUM ME- 
CHANICS: = 3 



In this appendix we explicitly demonstrate the relation between the quantum field the- 
oretic hamiltonian and the quantum mechanics many-body problem. First, we apply the 
periodic conditions and integrate by parts in the first part of hamiltonian ([1]) 

L L 

J,.ia.,^i.WM^--jM^W^. (Bl) 



For = 3, the A^-particle eigenstate | $) in ([7]) reads 

L 

|$) = j dXidX2dX3e~'^'^x{Xi,X2,X3)'i'^{Xi)'^^{x2)'i'^{x3)\0), Xi<X2<X3. (B2) 


L 

We consider operations H on this state. Let A = J dx[dx'^\x)][dx'i! (x)] \ $), then 



L 

I 



L 



A=-J dxdxidx2dx3e 'a^^xla^i, a;2, X3)(9^*'^(x))^(x)^^(xi)^^(x2)^^(x3) 1 0) 



L 

^dxdxidx2dx3e-'''i^xixi,X2,X3)idl'^\x)) (e-'''"'(^'^i)^^(xi)$(x) ¥{X2)'^\X3)\0) 

L 

' dxdxidx2dx3e~''^'^x{.xi,X2,X3){dl^\x))5{x - xi)'^\x2)'^\x3)\Q) 





L 





L 



-/..........3e-f"X-(......3)(a^*'(.))e-'-'*'(.0*(.)*'(..)*'(.3)|0), (B3) 



Here we succeeded in exchanging the position for and by means of the commu- 

tation relations ([2]). In flB3l) . the first term is the right term. We continue this procedure 
for the second term until \I'(x) is moved to the right hand side of \E'^(x3). Thus 

L 

A= -jdxdxidx2dx3e~''^'^x{xi,X2,X3){dl^\x))5{x - Xi)\P^(x2)^^(x3) |0) 



L 

-fdxdxidx2dx3e'^^^^^^'"^^'^^^^xLxi,X2,x^{dl^\x)^^^^ 
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L 

dxdxidx2dx3e~'^'^x{xi, ^2, X3){dl'^'' {x))6{x - Xi)\&'^(x2)\l/'^(x3) |0) 




L 



L 

- dxdxidx2dxse~'^''x{xi,X2,Xs){dl'^''{x))6{x - a;i)^^(x2)^^(a;3) 1 0) 









L 




L 



-/(ia;dxida;2c/x3e-'(i'^+'^"'(^'^^));((xi,a;2,a;3)(a2^ta;F^(a;i)5(a;-X2)^ta;3)|0) 







-Jdxdxidx2dxse-'^'^^+'''^^''^^'^+^'^^''^''^^^xixi,X2^^ 



The results ^(x)^^(x3) = e-'''"'(^'^3)^t(a;3)^(a;) + ^(a; _ 3,3) and ^(x) | 0) = are used in 
the last step. Using integration by parts and the integral property of the 5-function, we then 
have 

L 

A = -jdxidx2dx3e-'^''{dl^x{xi,X2,X3))'^\xi)¥{x2)'i!\x3)\0) 


L 

-/dxidx2rfx3e-^(t«+'^"'(^-2'^i»(9^^;^Xi,X2,X3))^^(x2)^t(xi)*t(a;3)|o) 





L 







dxidx2dx3e-'^'^^+^'"^''''^'^-^^'"^''^'''''^\dl^}ixi,x^^^^^ |0). 



dxidx2dx3e-''2-{{dl+dl+dl)x{xuX2,X3))^^Xi)^^{x2)^^x,)\0). (B4) 







Here the commutation relations ([2]) are used to form the pattern \l/"'"(xi)\l/^(x2)\l'^(x3) by 
moving around \E'"''(x2) and '^''{x^). 

We now turn to the second part of the hamiltonian (Q. Accordingly, let B = 
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L 

J dx'^^{x)'^''{x)'i{x)'i{x) I $), then 





L 

B 





= J dxdxidx2dx^e ^2, X3)*'''(x)*^(x)*(x)*(x)*'''(xi)*^(x2)*^(x3) | 0) 


L 

dxdx-rdx2dx3e-'^^x{x^X2,X3X^\x)y^{x){^\xi)^{x)e-'^'^^'^^^^^ 



L 

= / dxdXidX2dX3e-''^'^^^^'^^''^'''^\{xi,X2,X3){<i/\x) f^{x)¥{Xi)^ 





L 



+ J dxdxidx2dx3e '^"xi^i, ^2, X3){'if''{x)f'if{x)5{x - a;i)*^(x2)*^(a;3) 1 0) 



L 

=^Jdxdxidx2dx3e-'^'^^+^'^^'''''''^^x{xi,X2,X3X'i/\x)f^{x)^^ 

dxdxidx2dx3e~^^''x{xi,X2,X3){'^\xi)f'^{xi)'^\x2)'^\x3)\0) 





L 





L 

+Jdxdxidx2dx3e'^'^^+^'^^''''''^^x{xi,X2,X3X'if\x)f'i/{x)'i/\xi^^^ 



L 

+ [ dxdxidx2dx3e-''^''x{xi,X2, X3){'i/^{xi)y (*^(a;2)*(a;i)e-^'^'"(^i'^^) + 5{xi - X2)) ^^(2:3) 1 0) 





L 



^Jdxdxidx2dx3e'^'^''+'''^^''''''^^+'''^^''^'''^x{xi,X2,^^^ - X3) \ 0) 



L 

+ y"da;idx2dx3e-^(t''+'^"'(^2'^i^)x(2;i,2;2,a;3)(*^(a;2))^*(a;2)*^(a;i)*t2^3) | 0) 


L 

+ [ (ia;(ia;i(ia;2rfa;3e-'(i'^+'^"'(^i'^2))x(^i,^2,a;3)(*n^i))^*n^2)*(a:i)*H^3) |0) 





L 



+ J dxdxidx2dx3e '^'^x{xi,X2,X3){'if\xi)y5{xi - X2)'if\x3)\0) 
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L 

+ /"(ia;i(ia;2(ia;3e-'(5«+«"'(^2'^i))x(a;i,a;2,a;3)(^^(a;2))^^(a;2)^^(a;i)^ta;3) I 0) 




L 



+y"dxidx2dx3e-^(i'^+'^"'(^i'^2»x(2;i, X3)(*t2^i))^*t2^2)(*t2;3)*(a;i)e"'''"'^''^'''^^+5(xi-X3))| 



L 

+ y"cixicix2cix3e-'i'^x(2^i,a:2,a;3)(*^(a;i))^(^(a;i -2;2)*^(a;3)|(|. (B5) 



Denote these four last terms by 64, 63, 62 and h\ in order from top to bottom. 

We now discuss these terms case by case. The simplest case is 61. With the help of the 
property of the b function it becomes 

L 

61 = y" dxxdxidxj,Q-''^^''x{x\^X2, x^)?)[xx - X2)^\xx)^\x2)^\x^) I 0) (B6) 


as expected. By using \E'(a;) | 0) = 0, we have 

L 

62 = y"rfa;ida;2cix3e-'(i'^+'^"'(^i'^^))x(a:i,a:3a;3)5(a;i-a;3)(^\a;i))'^ta;2) |0) 


L 

= y"dxidx2dx3e-'(i«(^i'^^»x(^i, ^% x^)8{x^-x^)^\x^)¥{x^)^\x2) \ ^ 


L 

= /dxidx2(ix3e-^(t''+''^(^i'^^^)x(2^i,2:a2:3)(^(a:i-X3)*txi)e''^'"(^^'^2)*V2)*t2;3)|Q 





L 



y"(ia;irfa;2(ia;3e-^t'^x(a^i,a;2,a;3)(5(a;i-a;3)^'^(a;i)^^(x2)*\x3)|0). (B7) 







The term 63 is more complicated because we have to move '^{x'l) to the right hand side of 
'^'^{x-i)^ namely 

L 

dxicix2C^X3e-^''(i+"'(^2'^i))x(^i,2:2,a;3)(*^(2;2))^ (*^(2;i)*(x2)e-''^^(^2'^i) + -X2)) ^^^3) I 0) 

Ux^dX2dXz^'^^''"''^'''"^''^^'''^\{xx,X2,X^\'^\x2)^^^ 
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L 









L 



+jdxidx2dx^e'i^x{xi,X2,x^)5{xi-X2)^\xi)^!\x2)^\x^)\Q) 


L 

=jdxidx2dx^e-'^'^''+^'''"^''^^''^^\{xi,X2,x^)^Kx2)^Kx-^^^ |0) 



L 

+jdxidx2dx^e'^i^x{xi,X2,x^)5{xi-X2)^\xi)^!\x2)^\x^)\Q). 



Now using commutation relations ([2]), we move to the left side of \I'^(a;2) such that 

L 

b3= dxidx2dx3e~'^'^x{xi,X2,X3)S{x2 - X3)^'^(Xi)^^(x2)^'''(x3) |0) 





L 



+ dxidx2dx3e '^''x{xi,X2,X3)S{xi-X2)'^\xi)'$\x2)'^\x3) \ 0) . (B8) 





Now consider the term 64. Similar to the term 63, we have to move "^{xs) to the right hand 
side of \l/^(x2) in order to form the state \l/"'"(xi)\l/'^(x2)\l/^(x3) | 0). We thus have 

L 

bi=[dxidx2dx3e~'^^'^+'"^''''^''''^^'"^^''''^''^^xi^^ 





L 





L 





L 



rfxirfx2rfx3e-'"(2+2-(-3,xi)+-{x3,x2))^(^^^3,2^^3)(^t(a;3))^t(a;^)^(^^ 

+yc/xic/x2t^X3e-'"(i+"'("^'"^)+'"("^'"2))x(a^i,2:2,2:3)'5(x3-xi)(^tx3))'^tx2)|0) 


dXidX2dX3e-'''^'2+2'^^^^'^'^+^^^''^'^^^^{Xi,X2,X^^^ 
L 

+ /dxidx2rfx3e-''^(i+"'(^^'^^))x(a;i,a;2,X3)5(a;3-xi)^ta;i)*^(x3)^tx2)|0) 
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L 







L 





L 





L 





L 

+y"rfa;idx2rfx3e-^t'"x(a;i,a;2,a;3)5(a;3-a;i)^ta;i)^ta;2)*ta;3)|0) 



Therefore 



L 



+yrfxirfx2C?X3e"^2'"x(a;i,X2,X3)5(x3-Xi)^ta;i)^ta;2)^ta;3)|0). (B9) 


Summing the individual resuhs for bi flB6p . 62 flB?^ 63 fIBSP and 64 flB9p gives the result 



L 



i? = 2 y dXidx2dx3e '2«;i^(xi, X2, X3)(5(Xi-X2)+5(Xi-X3)+5(x2-X3))\E''^(Xi)\l'"'"(x2)^'^(x3) | 0). 

(BIO) 

Finally, we have 



1$) = ^A+\g,oB 
2m 2 



3 o9 3 



y dxidx2dx3e 'i'^ f 1^ + -^(^^ - 







j=l l<i<j<3 



X^t(;Xi)^^(x2)^"^(x3)| 0) 

E I $) (Bll) 



provided x satisfies 



i2 3 rj2^ 3 



-^J2d^ + 3iD ^^^^ - ^^■)^ = (^12) 

i=l « l<«<i<3 

which is the desired three-particle quantum mechanics problem. Here x = x{^i}X2,X3) and 
E is the eigenvalue. 
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APPENDIX C: PERIODIC BOUNDARY CONDITIONS 



Application of the periodic boundary condition on the Bethe wavefunctions for = 3 
gives 

X{xi = 0, X2, X3) = e-it("'(^i'^2)+«.(xi,x3)+«,(x2,x3)) J2 A{kpi, kp2, fcp3)e'(*=^2x='+fcp3^3) (qI) 

p 

X{x2,x3,xi =L)= e-it«^2'^«)+"'(^2'^i)+"'(^3'^i))^A(fcQi,fcQ2,fcQ3)e^(^'«i^^+^^^^^ (C2) 

Q 

We thus have 

,if (2«;(zi,X2)+2«;(xi,X3)) ^ ^(^^^^^ ^^^^ ^^^^ )gi(fcp2^2+fep3^3) 



6 2' 



L) 



(C3) 



Namely 



p 

Q 



(C4) 



according the order xi < X2 ■ ■ ■ < xn- The condition (IC4p gives six equations. Only three 
of them 



A(A;2A;3A:i) 



J2k 



J2k 



A{kik3k2) 



{ki 


-k2 


+ ic'){ki 


-k3 


+ ic') 


{ki 


-k2 


- ic'){ki 


-ks 


-ic') 


{k2 


-k. 


+ ic'){k2 


-ks 


+ ic') 


ik2 


-k. 


- ic')(^2 


-ks 


-ic') 


{ks 


-h 


+ ic')(A;3 


- k2 


+ ic') 


{kz 


-k. 


- ic')(fc3 


- k2 


-ic') 



(C5) 
(C6) 
(C7) 



A{k,k2k3) 

are independent. Similarly consideration of the case x{xi^X2, X3) = e^'^'^xi^i^ X2 + L, X3) gives 
the same set of equations. 
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